Reversible enhancement of the magnetism of ultrathin Co films by H adsorption 
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By means of ab initio calculations, we have investigated the effect of H adsorption on the struc- 
tural, electronic and magnetic properties of ultrathin Co films on Ru(OOOl). Our calculations predict 
that H occupies hollow sites preserving the two-dimensional 3-fold symmetry. The formation of a 
complete H overlayer leads to a very stable surface with strong H-Co bonds. H tends to suppress 
surface features, in particular, the enhancement of the magnetic moments of the bare film. The 
H[-induced effects are mostly confined to the Co atoms bonded to H, independent of the H coverage 
or of the thickness and the structure of the Co film. However, for partial H coverages a significant 
increase occurs in the magnetic moment for the surface Co atoms not bonded to H, leading to a net 
enhancement of surface magnetism. 

PACS numbers: 68.43.Bc,73.20.Hb,75.70.Ak,75.70.Rf 



I. INTRODUCTION 

Hydrogen is an ubiquitous element, present even in 
ultra-high-vacuum chambers, although very difficult to 
detect. It shows a large sticking probability in most tran- 
sition metal surfaceai, among them Co. Adsorption of H 
on transition metals was intensively studied in the past 
decades, mainly motivated by the understanding of the 
action of residual gases in catalytic and embrittlcmcnt 
processes. Much less attention has been paid to the influ- 
ence of H on the magnetic properties, though it is known 
to have a remarkable effect. For example, the adsorption 
of H2 molecules on the surface of magnetic transition 
metals produces a reduction of surface magnetism^. In 
addition, H coverage shifts the critical thickness for the 
spin reorientation transition (SRT) of Ni films^, and a re- 
versible manipulation of the magnetic exchange coupling 
of Fe/V superlattices can be achieved upon loading the 
V layers with H^. 

Previous experimental results characterizing the struc- 
ture and energetics of H adsorption on different Co 
surfaces and nanostructures with two-dimensional (2D) 
3-fold symmetry, as Co(OOOl)^ or Co nanoislands on 
Cu(lll)§., indicate that H prefers adsorption sites pre- 
serving the symmetry of the 2D lattice. This is in agree- 
ment with ah initio calculations of the hep Co(OOOl) 
termination^. Even for Co(lOlO) the formation of super- 
structures leading to adsorption sites of 3-fold symmetry 
has been reported^. The evaporation of molecular H2 on 
different Co surfaces results in the dissociative chemisorp- 
tion of H atoms up to a monolayer (ML) completion^. For 
higher doses, molecular II2 physisorption occurs. The 
desorption of H from Co nanoislands on Cu(lll) takes 
place locally and leaves the Co surface as it was prior 
the H evaporation, in contrast to the substrate exposed 
areas^. This surface character of the H-Co interaction is 
also supported by total energy calculations, which show 



a negligible weight of the H induced states even in sub- 
surface Co atoms^i^. In addition, no dependence of the 
binding energy on the H coverage is obtained^. More- 
over, H subsurface sites are highly unstable, indicating a 
low probability for H diffusion to the bulk. 

As compared to other transition metal surfaces, there 
are only a few studies of the influence of H on the elec- 
tronic and magnetic properties of Co surfaces and thin 
films. Spectroscopic measurements indicate a quench- 
ing of the sp-like surface states and a shift of the d sur- 
face states downward in energy upon H adsorptiortiSiii. 
Moreover, the Co electronic properties can be tuned in- 
dividually controlling the amount of adsorbed Hi^. The 
calculated density of states (DOS) also indicates that H 
states appear below the d valence band (VB)S'. When 
spin-polarization is included, it can be shown that this 
results in the quenching of the surface-induced enhance- 
ment of the Co spin moment, at least for 1 ML of H 
coverage^. 

Recently, we studied the magnetic properties of ultra- 
thin Co films grown on Ru(0001)i^. This is a unique 
system showing a double SRT linked to the completion 
of Co atomic layers. The SRT is intimately related to 
the structural distortion of the ultrathin Co film. Both 
Co and Ru are hep metals, but with very different lattice 
parameters: 2.51 A and 2.71 A, respectively. While 1 
ML of Co grows epitaxially, the second layer starts com- 
pressing towards the Co lattice, leading to an intermedi- 
ate 2D lattice parameter of 2.6 A. In this work we study 
the influence of H adsorption on the structural, electronic 
and magnetic properties of these ultrathin Co/Ru(0001) 
films, with special emphasis on the 2 ML thick ones which 
show perpendicular magnetic anisotropy. The study of 
the magnetic anisotropy will be the subject of a forth- 
coming paper. Our results indicate that H bonds strongly 
to the surface and tends to form a complete overlayer 
leading to a very stable surface system. Partial H cover- 
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age, although thermodinamically stable, presents slightly 
smaller energies of adsorption. As we will show, while 1 
ML of H quenches the surface magnetism, lower H cov- 
erage may enhance it. This can be a crucial factor re- 
garding the role of residual H in the measured magnetic 
properties of Co surfaces. 



II. COMPUTATIONAL METHODS 

We have performed ab initio calculations within the lo- 
cal density approximation (LDA), combining two differ- 
ent approaches. First, we have performed an exhaustive 
search of the equilibrium adsorption positions of H using 
slab models, a plane-wave basis set and the projector- 
augmented wave (PAW) methodic as implemented in 
the Vienna Ab-initio Simulation Package (VASP)i^. For 
the most stable configurations, we have then determined 
the electronic and magnetic properties within the fully- 
relativistic framework of the Screened Korringa-Kohn- 
Rostoker (SKKR) methodic, which allows for a smooth 
matching of the surface region to the semiinfinite bulk 
and lacks of the quantum size effects inherent to slab 
models. 

For the calculations with VASP, an energy cutoff of 350 
eV has been used. The structures arc modeled by peri- 
odically repeated slabs of 6 to 10 Ru layers covered by 2 
Co MLs on one side, and adding 1/4 ML or 1 ML of H on 
top of or underneath the outermost Co layer. A vacuum 
region of 12 A is left between both slab surfaces which 
prevents any spurious interaction between them. This 
model also guarantees the recovery of the bulk proper- 
ties at the inner Ru layers of the slab. The search for the 
most stable H adsorption sites has been performed start- 
ing from different H geometries and allowing full relax- 
ation of the atomic positions of H, Co and the outermost 
Ru layer. A uniform sampling of the k-mesh centered 
in r has been used, and convergence for the most stable 
structures has been obtained using k-samplings of up to 
(15 X 15 X 1) for (1 X 1) geometries and (8x8x1) for 
the (2 x 2) unit cells corresponding to adsorption of 0.25 
ML of H. Under these conditions, convergence in total 
energy differences below 1 meV is achieved. 

By using the optimized geometry with the most stable 
H adsorption sites, the detailed electronic structure and 
the magnetic moments were calculated in terms of the 
fully relativistic SKKR methodJ^ The SKKR formalism 
provides, in a natural way, layer-resolved physical quan- 
tities such as charges and magnetic moments. It, thus, 
manifests a perfect tool for the analysis of the H-Co in- 
teractions, and of the influence of the Co film thickness in 
the H induced properties. For this purpose, we have ex- 
tended our calculations to films of 3 to 6 Co MLs. Details 
of the SKKR calculations for the Co/Ru system can be 
found elsewherei^iii. The H atom is modeled by using an 
atomic sphere of radius 0.59 A, placed at the equilibrium 
distance to Co determined by VASP. 



III. STRUCTURAL PROPERTIES 

The stable adsorption site of H on 2Co/Ru(0001) has 
been explored starting from different geometries, de- 
picted in figure [T] top, bridge, hollow (either fee or hep) 
and off-symmetry positions, both at surface and subsur- 
face sites. In addition, for 1 ML coverage, the influence of 
structural modifications present in the Co film has been 
considered, such as the variation of the 2D lattice pa- 
rameter, or the existence of stacking faults which alter 
the Co/Ru stacking sequence from the hep AB/AB to 
BC/ABia. 




C) Co • H 

FIG. 1: (Color online) High symmetry sites for the adsorption 
of H on hep Co. 



A summary of the results is compiled in Table ID In all 
cases, the most stable adsorption site corresponds to H at 
fee hollow surface positions. The table shows the energy 
difference with respect to the equilibrium configuration 
only for the most stable structures. The rest of geome- 
tries considered, including on-top positions, provide total 
energies larger by more than 1 eV. Our results are in good 
agreement with previous ones relative to different Co sur- 
faces and thin films^i^: H tends to stick at the surface 
trying to preserve the 3-fold symmetry of the Co lattice, 
independent of the detailed Co film structure or the H 
coverage. The energy difference between adsorption sites 
provides an estimate of the diffusion barriers for H in the 
surface. Although the trends are similar for all surfaces 
studied, partial H coverage or alteration of the stacking 
sequence at the Co film may modify the diffusion barri- 
ers significantly. Regarding subsurface sites, the barrier 
for the diffusion of H into the Co film decreases as the H 
coverage increases, but it is always large enough to sup- 
port the low probability of the penetration of H into the 
film, in good agreement with the experiments of thermal 
desorption through formation of H2^. 

The results indicate that the adsorption of H is gov- 
erned by the local surface formation of strong H-Co 
bonds. Supporting this, for all structures considered the 
H-Co distance remains around 1.7 A, independent of the 
H coverage and the Co film structure. The largest vari- 
ations of the H-Co distance are found to depend on the 
H adsorption site, where H is coordinated to a different 
number of Co atoms: from 1.5 A for the on-top surface 
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adsorption (coordination 1), to 1.9 A for hollow subsur- 
face sites (coordination 6). 

TABLE I: Energy (in meV) per H atom with respect to the 
most stable fee hollow site for different H adsorption sites on 2 
ML Co on Ru(OOOl), both for 1 ML and 0.25 ML H coverage. 
The dependence on the in-plane lattice parameter a2D and 
the stacking sequence of the Co film are shown for the high 
H coverage. 



Co film LI adsorbate 

&2D (A) stacking hep bridge subsurf. 



1 ML 2.7 


AB 


27 


289 


> 420 




BC 


24 


708 


> 430 


2.6 


AB 


27 


278 


> 480 




BC 


108 


708 


> 550 


0.25 ML 2.6 


AB 


45 


236 


> 1000 



A further insight can be obtained comparing the ener- 
getics and structure of the equilibrium configurations of 
the Co film covered by different amounts of H. From now 
on we will restrict our considerations to the experimental 
value of a2D, 2.6 A. The adsorption energy per H atom 
is defined as: 

Eads = —^r;—[EHCoRu — EcoRu — NhEh] (1) 

where Nu refers to the number of adsorbed H atoms per 
unit cell, EhcoRu and EcoRu respectively to the total en- 
ergies of the slabs with and without adsorbed H and Eh 
to the energy of isolated H. The adsorption energies for 
the different H coverages are reported in Table [III Their 
positive sign indicates that adsorption is favorable in all 
cases. When comparing with the formation energy of the 
H2 molecule, 2.44 eV/H atom, we can see that both H 
coverages are highly stable against desorption, in good 
agreement with the experimental evidence and the pre- 
vious theoretical calculations for similar Co surfaces^i^. 
An alternative approach to determine the stability of the 
H covered surfaces comes from the evaluation of the work 
function, $, also provided in Table as obtained with 
the SKKR method. The increase of $ with H coverage, 
specially for the complete H overlayer, corroborates the 
large sticking coefficient of H at Co surfaces even for ul- 
trathin films, and is in good agreement with the results 
obtained for Co films on Cu(OOl)^. 

Table also provides the variation of the interlayer 
spacings in the Co film. As expected from the reduced 
atomic volume of Co compared to Ru, even for the bare 
surface there is a significant contraction of the Co-Co 
interlayer distance, dcci trying to compensate the in- 
plane expansion. Accordingly, the Co-Ru distance, dc_R, 
is only slightly reduced. The large volume compensation 
effect suppresses the oscillatory relaxation of interlayer 
distances typical of transition metal surfaces, although 
an oscillatory behavior can be identified in the layer- 
resolved electronic properties for larger Co thicknesses. 
Upon H adsorption, the progressive filling of free bonds 



TABLE II: Adsorption energy and work function for 
2Co/Ru(0001), either bare or covered by H, together with the 
variation of the Co-Co and Co-Ru interlayer spacings relative 
to the unrelaxed Ru bulk, Adcc and Adcjj, respectively. The 
two values of Adcc for 0.25 ML coverage correspond to Co 
atoms unbonded/bonded to H. 

H coverage F.ai, (eV) $ (eV) Adcc (A) Adcj; (A) 

Bare 6.03 -0.29 -0.06 

0.25 ML 3.34 6.09 -0.30/-0.27 -0.07 

1 ML 3.41 7.02 -0.22 -0.09 



attenuates the surface induced relaxations, but still the 
volume compensation remains. Thus, the resulting re- 
laxation pattern is different from that of other H covered 
films where Co retains its 2D lattice constant^. In spite 
of this difference, we can conclude that H always acts as 
an attenuator of surface effects. It is worth noticing the 
slight corrugation of 0.03 A at the Co surface layer for 
1/4 ML H coverage, driven by the Co atoms bonded and 
unbonded to H. Remarkably, the surface induced con- 
traction of dec is enhanced for those unbonded. 



IV. ELECTRONIC AND MAGNETIC 
PROPERTIES 

Next we discuss the layer-resolved electronic and mag- 
netic properties of the 2 Co/Ru(0001) film with and with- 
out H coverage. The atom projected DOS and magnetic 
moments are shown in Fig. [5] and in Table IIIIl respec- 
tively. The DOS of the bare film shows the narrowing in- 
duced by the reduced coordination at the surface plane. 
In the case of Co, with the majority spin band almost full, 
this causes an increase of the magnetic moment. As H 
adsorbs, the surface effects start to attenuate, and, upon 
completion of 1 ML of H, the magnetization of the two 
Co layers becomes almost uniform. The H-states localize 
at the bottom of the d valence band (VB), and show a 
significant spin polarization induced by the hybridization 
with Co, even though the net magnetic moment of H is 
almost zero. In turn, H-induced features can be identi- 
fied in a region dominated by sp states for the surface Co 
atoms bonded to H. The weight of these features for the 
rest of Co atoms, including those at the surface layer for 
0.25 ML of H (not shown in the figure), is negligible. 

The short range of the H-induced effects can also be 
inferred from the magnetic moments of the surface Co 
atoms bonded and not bonded to H for 0.25 ML cov- 
erage, see Table IIIIl Noticeably, the partial H coverage 
results in a significant increase of the magnetic moments 
of the Co atoms unbonded to H. This is evidenced both 
in the fully relativistic calculations performed with the 
SKKR method and in those performed with slab mod- 
els. Though the arbitrary assignment of the atomic radii 
affects the actual values of the magnetic moments for 
each atom, the relative differences between the structures 
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calculated under the same conditions provide a reliable 
physical picture. Furthermore, the total magnetic mo- 
ment per slab, which is not subject to the assignment of 
the individual atomic radii, also supports the enhance- 
ment of the surface magnetism for partial H coverages: 
9.64 HB (bare), 10.88 fis (0.25 ML) and 9.72 (1 ML). 
This enhancement can be understood in terms of the 
strong local H-Co interaction, which limits the Co-Co in- 
teractions within the surface plane. Quite interestingly, 
similar trends for the orbital moments of the Co atoms 
upon H dosing can be inferred from Table IIIII From this 
we expect that the adsorption of H affects the magnetic 
anisotropy energy of the Co/Ru(0001) system in a non- 
trivial manner, as well. 

TABLE III: Layer resolved magnetic moments obtained with 
slab models (rrisiab), as well as spin (ms) and orbital (mz,) 
contributions to the magnetic moments calculated by the rel- 
ativistic SKKR method for the 2Co/Ru(0001) surface system 
covered by different amounts of H. The moments are given 
in /is and the labels of the atoms are the same as in Fig. 
[21 For 0.25 ML H coverage, the moments of the Co atoms 
unbonded/bonded to H are distinguished. 





Ru 


Co 


Cos 


H 






-0.05 


1.31 


1.69 




bare 




-0.07 


1.30/1.23 


1.76/1.54 


-0.02 


0.25 ML 




-0.02 


1.27 


1.27 


-0.01 


1 ML 


ms 


-0.03 


1.51 


1.75 




bare 




-0.01 


1.56/1.50 


1.94/1.75 


0.04 


0.25 ML 




-0.01 


1.52 


1.58 


0.01 


1 ML 




0.00 


0.08 


0.11 




bare 




0.00 


0.08 


0.13/0.10 


0.00 


0.25 ML 




0.00 


0.08 


0.07 


0.00 


1 ML 



Finally, we discuss the H induced effects for thicker Co 
films. To this end, we performed SKKR calculations for 
Co films of up to 6 MLs covered by 1 ML of H. Figures [3] 
and [4] show respectively the layer resolved values of the 
spin (m^) and orbital (m^) contributions to the mag- 
netic moment for the different Co thicknesses, together 
with those at the corresponding bare surfaces. Upon H 
adsorption there is always a strong reduction of both ms 
and m^ at the surface Co atoms. On average, the reduc- 
tion is around 0.25 fiB for ms and 0.04 for m^. By 
contrast, at the subsurface Co layer the reduction of the 
magnetic moment is low, becoming negligible for the lay- 
ers underneath. Although the specific values of ms and 
m^ at each layer involve in a complex way the oscillatory 



layerwise variation of the electronic properties, the effect 
of H does not depend on the Co thickness. Furthermore, 
it is very similar for both components of the magnetic 
moment. 

Not shown in the figure, in all cases the H atom ac- 
quires a ms of 0.01 /jb, and no orbital polarization. Sim- 
ilarly, the components of the moment induced at the in- 
terface Ru atom are around -0.02 /is (ms) and -0.001 
fiB (hil) for all Co thicknesses, without any significant 
variation as compared to the bare film. The layer re- 
solved DOS (not shown here) are consistent with this 
conclusion, with H induced features almost restricted to 
the surface Co layer in each case. Thus we can conclude 
that the H induced effects are essentially restricted to the 
surface and independent of the Co thickness. 



V. CONCLUSIONS 

In terms of ab initio calculations, we performed a 
detailed theoretical analysis of the stability, electronic 
structure and magnetism of Co/Ru(0001) films upon H 
adsorption. Our main observations are as follows: H ad- 
sorbs on Co at the hollow sites by preserving the threefold 
symmetry of the system. The adsorption is cndothermic 
for all H coverages, with the largest binding energy corre- 
sponding to the complete H overlayer. Furthermore, the 
probability of H desorption or segregation is low. Con- 
cerning the structural, electronic and magnetic proper- 
ties, the effect of H has a local surface character and 
is almost independent of the detailed Co thickness and 
structure. H always attenuates surface effects. However, 
for partial H coverages, this leads to the enhancement 
of the surface induced features for Co atoms not bonded 
to H, and in particular to an enhancement of the net 
surface magnetic moment. Our results point to the im- 
portance of residual H in magnetic measurements, and 
to the possibility for the reversible manipulation of the 
properties of ultrathin Co films and surfaces by tuning H 
adsorption. 
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FIG. 2: DOS projected on the H, surface Co (Cos), interface Co and Ru atoms of the 2Co/Ru(0001) surface system either bare 
(left) or covered by 1/4 ML (middle) and 1 ML (right) of H. The Cos DOS for 1/4 ML H coverage corresponds to Co atoms 
bonded to H. Energies are related to the Fermi level. 
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FIG. 3: (Color online) Layer-resolved Co spin magnetic mo- 
ment for Co films of 2 to 6 MLs (from top to bottom), either 
bare (empty circles) or covered by H (solid circles). Layers 
are numbered from the Ru interface. 




FIG. 4: (Color online) Same as figure[3]for the orbital contri- 
butions to the magnetic moment. 



